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Abstract 

We have measured the electrodynamic properties of Yi_ x Pr x Ba2Cu307 single crystal 
thin films as a function of temperature using coherent THz-time-domain spectroscopy. 
We obtain directly the complex conductivity a = o\ + io~2, the London penetration 
depth Al, the plasma frequency u p , and the quasiparticle scattering rate 1/r. We 
find that 1/r drops exponentially rapidly with T below the critical temperature in 
all the superconducting samples, implying that this behavior is a signature of high- 
T c superconductivity. The plasma frequency decreases with increasing Pr content, 
providing evidence that Pr depletes carriers, leaving the CuO planes underdosed. 
Both the conductivity in the THz region and the dc resistivity yield evidence for the 
opening of a spin gap above T c . 
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The large body of vital spectroscopic data on the high-T c superconducting compound 
YxE^CusC^-^ (YBCO) measured in the far infrared (FIR) as well as the GHz-region has revealed 
many interesting features. Among them are: the presence of one or more energy gaps in 
the superconducting state, a temperature dependent penetration depth that is indicative of d-wave 
pairing || , and a quasiparticle scattering rate which is proportional to frequency in the normal state 
|]T[ and which drops by four orders of magnitude at T c [|J. 

The new technique of THz-time-domain spectroscopy bridges the gap between the FIR and the 
GHz regime. THz measurements of YBCO display a peak in the temperature dependence of o\ below 
T c , which appears to be related to a temperature-dependent scattering rate 1/r rather than to BCS 
coherence effects |J. Other experiments using the same technique have revealed that the penetration 
depth Al, has a temperature-dependence inconsistent with BCS theory ||. 

Here we report results of a THz investigation into the Y!_ x Pr x Ba2Cu307 (YPrBCO) system. 
These measurements are of great interest because samples covering a wide range of critical tem- 
peratures can be obtained simply by varying the Y/Pr ratio 0. Also, because YPrBCO is fully 
stoichiometric in the oxygen content, it is more homogeneous than oxygen-depleted YBCO and does 
not change from an orthorhombic to a tetragonal lattice upon decreasing T c . Several mechanisms 
have been proposed for the T c suppression, including a magnetic pair breaking effect || and depletion 
of holes in the Cu0 2 planes ||. 

Recently, the opening of a spin gap at temperatures above T c in oxygen-deficient YBCO has been 
deduced from various experiments e.g. dc resistivity and Hall coefficient [|TU|,|TTf, neutron scattering 



and NMR measurements |I3|. However, until now only one electro- magnetic measurement on 



underdoped YBCO has been linked to the presence of a spin gap Q. In this paper, we report evidence 
for a spin gap obtained from measurements both at dc and at THz frequencies in YPrBCO, which is 
consistent with Pr depleting holes from the CuO planes ||. 

Our spectroscopic technique involves a coherent time-domain measurement of a ps microwave 
impulse E(t) transmitted through the sample ||14|| . A Fourier transform yields the complex transmis- 
sion spectrum t(u) and the complex conductivity o~(lj), without the use of Kramers- Kronig analysis. 
The London penetration depth A^(T), the plasma frequency uo p in the clean limit, and the scattering 
time t(T) are obtained with the aid of a two-fluid model. 



The microwave source is a biased 30 /im transmitting antenna fabricated on low temperature 
grown GaAs, triggered with ~100 fs optical pulses from a colliding-pulse modelocked dye laser. The 
emitted microwave pulses have spectral components spanning the 0.1-1.0 THz spectral region which 



is difficult to access with conventional electronics | I5| . The receiver is a 30 fim antenna, fabricated 



on ion-implanted silicon-on-sapphire and gated with a second pulse from the laser. The receiver 
photocurrent is proportional to the incident microwave field. 

We investigated five YPrBCO samples having Pr composition x = 0.0, 0.2, 0.3, 0.4, and 1.0. The 
films are grown by pulsed laser deposition |16[ onto NdGaC>3 substrates. The film thicknesses are 
approximately 150 nm (Table I). NdGaC>3 is the ideal substrate because it remains transparent and 
nondispersive over the entire spectral bandwidth of our pulses, as well as over the entire range of 
temperatures investigated here. The excellent lattice match of NdGaO"3 to YPrBCO is a prerequisite 
for a low density of misfit dislocations on the film-substrate interfaces. The substrates have a (001) 
orientation, yielding c oriented, twinned films. The critical temperatures T c dc (Table I) are determined 
from the temperature dependent resistivity, as measured with a four-point probe. 

To calculate the conductivity, we make use of the multiple reflection formula for the field trans- 
mitted through a layer of (complex) index 77.2, thickness d, bounded by media of index ri\ and n 3 : 

t 12 t 23 exp(in 2 (u/c)d) 
1 + r 12 r 23 exp{2tn 2 {uj / c)d) 

where tij = 2ni/(rii + rij) is the field transmission coefficient at the r/'th interface, and rjj = 
(rii — nj)/{rii + rij) is the field reflection coefficient. In our geometry, rii represents vacuum, 



n 2 — yl + ^ cr (^)/( e o^) is the index of the superconducting layer, and n 3 is the measured index 
of the substrate. As n 2 {uj /c)d <C 1, and n 2 ^> n 3 > 1 in our samples, Eq. ([[]) reduces to 

t{Uj) ~ l + n 3 + Z a{u)d (2) 

where Z is the impedance of free space. 

The effect of varying Pr content x (hereafter [Pr]) on the conductivity spectra <j(uj) at T = 50 K 
is similar to the effect of varying the temperature for a given alloy (Fig. |l|). The addition of Pr has 
at least two interrelated effects: a) The suppression of T c changes the partitioning between normal 
and superconducting carriers, b) The total number of carriers N (or their mobility) may be reduced. 



To the extent that the superconducting carriers make the largest contribution to 02, both factors a) 
and b) lead one to expect that at a given temperature, pure YBCO would have the largest a 2 , and 
it does. Samples with 20% and 30% Pr have smaller values of 02, since they are only slightly below 
their T c . The sample with 40% Pr has 02 ~ because it is above T c at 50 K. 

Only normal carriers contribute to o\ for uj 7^ 0, but now factors a) and b) compete. At 50 K, 
a\ decreases with [Pr], therefore the effect of a reduction in TV dominates the effect of the shift in T c 
which increases the fraction of normal carriers. The data for pure YBCO and 20% Pr have been fit 
to a frequency-dependent Drude form, and one can see that 1/r lies in our spectral range. For 30% 
and 40% Pr, we observe frequency independent o\ and can conclude that 1/r is large compared to 
1 THz. Pure PrBCO is a dielectric at 50 K, as seen by a conductivity proportional to frequency, i.e. 
a dielectric constant independent of frequency. 

Examining <t 2 at a fixed frequency, e.g. 480 GHz, we see that it is close to zero at high temperature, 
but rises sharply at the onset of superconductivity, thus providing an independent ac measurement 
of T c (Table I). 

In all of the superconducting alloys, below T c , <j\ displays a peak, which has been previously 
observed only in fully oxygenated YBCO using microwave techniques P,|17|| as well as measuring the 
thermal conductivity [18|]. It has been attributed to a sharp rise in the scattering time r of normal 



carriers below T c offsetting the decrease in the fraction of normal carriers, although coherence effects 
remain a possibility. For pure YBCO, the peak value is about 20 times higher than o"i(100 K). 
With increasing [Pr] the peak height in o~\ decreases but the temperature corresponding to the peak 
position does not shift significantly. 

The normal state behavior of our samples is particularly interesting because underdosed materials 
such as (124) YBCO and oxygen deprived (123) YBCO undergo a phase transition associated with 
the opening of a spin gap at a temperature To > T c . Evidence for the presence of a spin gap has 
been seen in neutron scattering |12|,|19j, as well as in dc resistivity measurements [|K|[n|]. Recent 



experimental |20| and theoretical || work confirms that YPrBCO alloys are also underdoped, i.e. 
superconductivity is suppressed, because holes are removed from the Cu02 planes. 

If the normal carriers couple strongly to spin fluctuations, the opening of a spin gap should be 
accompanied by an increase in the scattering time r, giving rise to an enhancement in <j\ below Td 



for u) < 1/t. For pure (optimally doped) YBCO, at 480 GHz, o\ shows only a single transition at 
T c (Fig. ||b). For the (underdoped) alloys, o\ has two transitions, one at T c , the other at a higher 
temperature which increases with [Pr]. To accentuate the two transitions, Fig. |2|b is shaded in the 
region bounded by T c , the experimental curve, and a dashed line representing l/(a + (3T) behavior. 

The higher transition temperature seen at 480 GHz matches that of a transition also observed in 
the dc resistivity (Fig. |3]a). The transition is manifested as a deviation from a linear T dependence. 
To show the deviation more clearly, the resistivity is first fit to a line p = a + (3T between 250 
and 300 K, and then the experimental values are normalized to the value determined by the line as 
follows: 

The normalized resistivity (Fig. |3|b) of the YPrBCO alloys clearly shows a transition occurring at a 
temperature above T c similar to that which has been observed in underdoped YBCO, and connected 



to the opening of a spin gap flQ|TT 



For the evaluation of Xl we are using a two fluid model of the form: 

€oLJ p T 1 ( xr \ i 



a(u) = - — x n -\ -y [-ir5(u) + — ) x s , (4) 

1 — iujt p X z L \ ujJ 

where uo p is the plasma frequency, and x n and x s are the fractions of normal and superconducting 
carriers, given by x n + x s = 1. We derive the penetration depth from a fit to Eq. (f|). In our frequency 
range, we can ignore the Drude contribution to 02 for temperatures which are more than a degree 
below T c . As x n = at T = in the two-fluid scenario, we have x s = (A_l(0)/ A^(T)) 2 . Shown in 
Fig. I is the measured (Al(0)/Ai / (T)) 2 vs. T/T^ c for all samples. Also plotted are theoretical curves 
from weak-coupling BCS theory [^l|] and the functional form 



(A L (0)/A L (T)) 2 = 1 - (T/T c r. (5) 
Power-law behavior in the dependence of A^ on T for T — > is one signature of "unconventional 



superconductivity", i.e. nodes in the energy gap in fc-space ||22|| . Our data are intermediate between 
BCS theory and a = 2, which is the exponent predicted for <i-wave pairing with impurity scattering 
23| . Microwave measurements at lower frequencies have indicated both a = 1 (in the T — >• limit) 



[PljMfl and a = 2 0. The Gorter-Casimir two- fluid model (a = 4) [p5|1 , as well as c?-wave pairing 
(a = 1) are clearly inconsistent with our data. 

Al(0) is obtained by extrapolating \l(T) to T = using a quadratic regression (Table I). Our 
value for pure YBCO, 168 ± 2 nm, is close to the 145 nm value typical for good samples [f|. Ax(0) 
increases with increasing [Pr] as the samples become less superconducting. 

The plasma frequency uo p (Table 1) is determined with the relation uo p = c/A^(0), which is valid 
in the clean limit, provided all carriers are superconducting at T = 0. The clean limit assumption 
becomes less valid upon substitution of Pr, therefore the values for the alloys are only estimates. 
Our value of 9500 cm -1 in pure YBCO is close to the 12000 cm -1 value indicative of good sample 
quality f|. 

The quasiparticle scattering rate is calculated by fitting Eq. (§) to the data for <Ji and using, for 
x n , the relation x n = 1 — x s = 1 — (Al(0)/A^(T)) 2 . We find that the scattering rate 1/r decreases 
slowly with decreasing temperature until T f=a T c , after which it decreases exponentially rapidly 
(Fig. ^|), possibly reflecting the opening up of a gap in the fluctuation spectrum. Similar behavior 
has been reported in YBCO samples |4j]. The scattering rate 1/r reveals an exponential decay below 
T c by more than two orders of magnitude for all samples with different [Pr] |26| . Thus we are led to 
conclude that this behavior is a universal feature of high-T c superconductivity. 

For T > T c , we find the surprising result that 1/r decreases with [Pr] implying that alloy scattering 
plays a negligible role, and that depletion of carriers in the Cu02 planes dominates the scattering. It 
should be noted, however, that the drop in 1/r with [Pr] may be due to an underestimate of ui p , if the 
clean limit assumption breaks down for the alloys. A more precise determination of the scattering 
rate dependence on [Pr] above T c is beyond the scope of this letter. 

In conclusion, we have presented complex conductivity experiments in the THz range on thin 
films of Yi_ x Pr x Ba2Cu307. <Ji{T) reveals the anomalous coherence peak for all superconducting 
alloys. Both the conductivity in the THz region and dc resistivity measurements provide evidence 
for the opening of a spin gap in the excitation spectrum of the underdoped samples. a 2 (u;,T) yields 
directly the penetration depth. For the superconducting YPrBCO alloys the plasma frequency u p in 
the clean limit decreases with [Pr] above T c due to a reduced number of carriers. The temperature 
dependent quasiparticle scattering rate l/r(T) shows an exponential drop below T c in all the alloys. 

6 



These observations fit into the picture of Pr suppressing superconductivity by depletion of the carrier 
concentration in the superconducting Cu0 2 planes. 
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FIGURES 

FIG. 1. a) d\{y) and b) oiiy) for Yi_ x Pr x Ba2Cu307 at T = 50 K, x = (squares), 0.2 (circles), 0.3 
(dotted curve), 0.4 (dash-dotted curve), and 1.0 (short-dashed curve). The solid and dashed upper curves 
are fits to a frequency-dependent Drude form (see text). 

FIG. 2. a) o"2(T) and b) 0"i(T) for the four Pr concentrations x = (squares), 0.2 (circles), 0.3 (solid 
triangles), 0.4 (open triangles) at 480 GHz. The dashed curves in b) show 1/p*. The vertical lines indicate 
T° c . For clarity, the data for x = 0.3 in b) has been multiplied by 1.5. 

FIG. 3. a) The dc resistivity p(T) and b) the normalized dc resistivity p*(T) for x= 0, 0.2, 0.3, and 0.4 
(solid, dashed, dotted, and dash-dotted curves). 

FIG. 4. Normalized London penetration depth Xl(T /T^ c ) for x = 0.0 (squares), 0.2 (circles), 0.3 (solid 
triangles), 0.4 (open triangles). Also shown are curves corresponding to the predictions of BCS theory (solid 
line), as well as to the functional form (|5|) for a = 1 (dash-dotted line), 2 (dashed line), and 4 (dotted line). 

FIG. 5. Quasiparticle scattering rate l/r(T) for x = (squares), 0.2 (circles), and 0.4 (triangles) at 
480 GHz. 



Table I 



X 


d (nm) 


Al(0) (nm) 


Up (cm x ) 


T c dc (K) 
(4-point probe) 


7? (K) 
(microwave) 


0.0 


155 


168 ± 2 


9500 


93 


92 


0.2 


134 


349 ± 10 


4560 


68 


72 


0.3 


170 


375 ± 7 


4244 


53 


59 


0.4 


170 


590 ± 15 


2693 


40 


41 
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